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ABSTRACT 

HESS J0632+057 is a variable, point-like source of Very High Energy (> 100 
GeV) gamma-rays located in the Galactic plane. It is positionally coincident with 
a Be star, it is a variable radio and X-ray source, has a hard X-ray spectrum, and 
has low radio flux. These properties suggest that the object may be a member 
of the rare class of TeV/X-ray binary systems. The definitive confirmation of 
this would be the detection of a periodic orbital modulation of the flux at any 
wavelength. We have obtained Swift X-ray telescope observations of the source 
from MJD 54857 to 55647 (Jan. 2009 - Mar. 2011) to test the hypothesis that 
HESS J0632+057 is an X-ray/TeV binary. We show that these data exhibit flux 
modulation with a period of 321 ± 5 days and we evaluate the significance of 
this period by calculating the null hypothesis probability, allowing for stochastic 
flaring. This periodicity establishes the binary nature of HESS J0632+057. 

Subject headings: gamma-rays: observations — X-rays: binaries — X-rays: indi- 
vidual(HESS J0632+057) 



Introduction 



There are cu rrently three confirmed T eV gamma-ray bina ries; PSRB1259-63 (lAharonian et al. 



2005at ). LS 5039 f kharonian et al. 12005 bL and LS 1+61° 303 fUlbert et al. Ibood : lAcciari et al 
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2008, 2009a ) as wel l as a fourth new GeV gamma-ray binary recently detected by Fermi 
( ICorbet et al. 1 1201 lh that could have associated TeV emission. Evidence at the 4.1a level 
for Very High Energy (VHE) emiss i on fro m the stellar mass black hole candidate Cyg X-l 
has been reported by I Albert et al. I (120071 ) during a single short flaring episode. Contained 
within each of these sources, a young and massive O or B star is orbited by a compact object, 
either a neutron star or a black hole. The VHE emission in TeV binary systems is powered 
either by pulsar winds driving shock acceleration or by mass accretion onto the compact ob- 
ject driving a microquasar jet. The compact physical size of the emission region leads TeV 
binaries to be one of the few classes of Galactic VHE sources which appear point-like to TeV 
instruments. Following the discovery of several TeV objects with no identified counterparts, 
HESS J0 632+057 was suggested to be a possible mem ber of this exotic class of TeV/X-ray 
binaries ( lAharonian et al. 1 120071 ; iHinton et al. 1 120091 ) . We are testing this hypothesis by 
using X-ray data to search for periodicity from orbital flux modulation. 



HESS J0632+057 was first detected as part of the HESS Galactic Plane Survey ( lAharonian et al 



20071 ). Located on the edge of the Rosette Nebula within the Monoceros Loop region at R.A. 



06h 32' 58.3", Dec. +05° 48' 20" (± 28" stat., 20" sys.), the spatial distrib ution of the source 



is con sistent with that of a point source (rms size < 2' at 95% confidence) (lAharonian et al. 



20071 ). The source is located on the edge of the 99% contour of the EGRET unidentified 
object 3EG J0634+0521, but no object at the position of HESS J0632+057 is listed in the 
Fermi LAT point source catalog. MWC 148, a BOpe star, is positionally coincident with the 
centroid of the HESS position. 



Hinton et al. I ( 120091 ) observed the region surrounding HESS J0632+057 on 2007 Septem- 
ber 17 with XMM, obtaining 26 ks of good data. The observation resulted in the identification 
of point source XMMU J063259.3+054801, which is positionally coincident with both MWC 
148 and HESS J0632+057, and within the 99% error circle of 1RXS J063258.3+054857. The 
source showed significant variability during this short observation (« 7.2 hours). During this 
observation the source showed a mean deabsorbed flux of (5.3 ± 0.4) xl0~ 13 erg cm -2 s" 1 
in the 1-10 keV range. The spectrum was well fit by an absorbed power-law model with a 
photon index of 1.26± 0.04 and a column density of Njj = 3.1 ± 0.3 x 10 21 cm -2 . 



Beginning 2009 January 26, [Falcone et al. I (j2010h initiated a campaign to obser ve 
XMMU J063259.3+054801 with the Swift X-ray telescope (XRT; burrows et al~l j2005[ )). 
reporting significant variabilit y in 0.3-10 keV flux. While flux variability was measured on 
day to month-long timescales, iFalcone et al. I (120101 ) found no evidence for periodicity and 
concluded that if periodicity was present, it was likely with P > 54 days. Additionally, the 
X-ray source has shown no e yidence of short tim escale periodic or quasi-periodic variability 
in the range 0.005 - 800 Hz faea fc Torres lb 111 ). 
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During December 2006 - January 2007 and December 2008 - January 2009 lAcciari et al. 



(I2009bf ) performed follow-up observations of HESS J0632+057 above 1 TeV with VERITAS, 
yielding flux upper limits well below the values published by HESS. Together, the VERITAS 
and HESS observations provide evidence for variability of the gamma-ray flux on time scales 
of months. More recently, both VERITAS and MAG IC have detected elevated TeV gamma- 



ray emission (lOng et al. 



flux reported by iFalcone et al. 



2011 



Mariotti et al. 11201 ll ) during the time period of elevated X-ray 



20111 ) in February 2011. 



W e find no report of optical flux modulation in the literature although lAragona et al. 
(120101 ) report temporal variations in the Ha emission line profile. HESS J0632+057 has 
shown a radio spectral index of a r = 0.6 ±0.2 and significant variabilit y at 5 GHz on month- 



long time scales around a mean flux of 0.3 mJy (ISkilton et al. 1120091 ). however no periodic 
variability could be detected in these data. This radio flux is much lower than the typical 
radio flux expected from a TeV blazar, making such a potential interpretation improbable. 
Recently, iMoldon et al. I ( 120111 ) announced the 1.6 GHz detection of a milliarcsecond scale 
source coincident with MWC 148. 

The point-like nature of the detected TeV source, the e xcellent positiona l coinc idence 



with MWC 148 (chance coincidence of ~ 10 4 according to lAharonian et al. I (120071 )). the 



location on the Galactic plane with a low radio flux, the X-ray binary-like spectral index, 
and the variable X-ray and gamma-ray emission are all facts that argue in favor of an X-ray 
binary in association with MWC 148. Confirmation of HESS J0632+057/MWC 148 as a 
TeV binary would add a new member to this short list of objects; whereas, refutation of 
the binary hypothesis would establish an equally interesting class of TeV source. The most 
direct way to test this binary hypothesis is to search for periodic emission signatures. In this 
paper, we report on recent monitoring data taken with Swift-XKT. 



2. The Observations 



The X-ray Telescope (XRT; iBurrows et al. I ( 120051 )) on the observatory (IGehrels et al. 



20041 ) was used to obtain sensitive observations in the 0.3-10 keV energy band. The typical 
observation duration was ~ 5 ks with ~ 1 week spacing between most observations. From 
16 May 2009 - 15 August 2009 and 10 May 2010 - 21 August 2010 the source was unob- 
servable due to its close proximity to the Sun. The total Swift-XKT data set includes 463 
ks of observations, extending from 2009 January 26 (MJD 54857.1) to 2011 March 27 (MJD 
55647.6), representing a baseline of T = 790.5 days. All of the observations were obtained in 
photon counting (PC) mode. All observations beyond MJD 54965 are being reported here 
for the first time. 
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3. Analysis 



We used the most recent versions of the standard Swift tools and the most recent cali- 
bration files available at the time of data processing. In particular, we utilized Swift Software 
version 3.7, FTOOL S version 6.10, and XSPEC version 12.6.0q. Light curves were gener- 
ated using xrtgrblc version 1.5. Source and background regions, bad CCD detector column 
correction, a nd point spread funct ion corrections, and filtering of data was all performed as 
described in iFalcone et al. I (120101 ) . Since these observations always resulted in XRT count 
rates ~0. 01-0. 08 c/s, there was no significant pile-up. 

Spec tral fitting was done with an absorbed power-law, with th e N H set to the valu e 
found bv lliinton et al. I ( 120091 ) with XMM. Following the technique of Falcone et al. I (1201 (J ), 
we calculated rate-to-flux conversion factors for groups of observations, defining the groups 
such that they contained enough counts to enable good spectral fits ( 20 counts/spectral bin) 
and isolated a single flux state of the source. Assuming that the observed energy spectrum 
does not change significantly while the source is in a particular flux state, the rate-to-flux 
conversion will apply for all observations in each group. The time-dependent flux-rate ratio 
was then applied to the rate light curve, shown in Figure 1. Error bars represent the count- 
rate error scaled by the rate-to-flux conversion factor. 



4. Results 

A spectral fit of the combined data set finds that the data is well fit by an absorbed 
power-law spectrum with column density set to the XMM value (Nh = 3.1±0.3 x 10 21 cm -2 ). 
This results in an average spectral index of 1.58 ± 0.06 and x 2 — 5 7.9 with 39 deg r ees o f 



freedom. This time averaged spectrum is consis tent with the result s of IFalcone et al. I (120101 ) 



and other TeV binaries such as LS 1+61° 303 (jSmith et al. 1 120091 ) . The X-ray light curve, 
constructed from the entire data set, is shown in Figure [TJ Given the observations, we define 
five different feature types in the variability; low level quiescent emission of ~ 10~ 12 erg cm -2 
s _1 , large emission peaks, which are 4-5 times brighter than the quiescent emission, small 
precursor peaks, which are approximately 2 times brighter than the quiescent emission and 
precede each large peak, small mid-phase peaks, which are approximately 2.5 times brighter 
than quiescent emission and appear approximately half way between the large peaks, and 
dips, which succeed each large peak and represent a decrease in emission by a factor of 
~ 1/3 below quiescent. By eye, these emission features appear to be spaced at regular 
intervals, leading us to believe that their positions reveal an underlying periodicity in the 
source. To confirm the presence of periodic modulation, we apply an autocorrelation analysis. 
The period of the modulation is determined with a peak fitting method, and its statistical 
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significance is established using MC light curve simulations. 



4.1. Autocorrelation 



An autocorrelation analysis is often used as a test for repeating patterns i n astronomical 
time series, as it requ ires no prior assumptions about the lightcurve profile (jCiprini et al. 
20071 ; Kie et al. 1120081 ). The Discrete Auto-Correlation Function (DACF), in particular, can 
be used to study the level of auto-correlation in unevenly sampled data sets without any 
interpolation or addition of artificial data points (Ed elson &: Krolik Ill988l ). Figure [2] shows 
the result of applying the ^-transformed DACF (I Alexander 1 1 19971 ) to the Swifl-XBT dataset. 
We show both the autocorrelation value, and the autocorrelation divided by its error, which 
provides an estimate of the significance of the autocorrelation measurement. A clear peak is 
visible at a time lag of ~ 320 days, as well as the expected peak at twice this period. 



4.2. Peak Fitting 

Peak fitting is a simple and robust method for measuring periodicities in data, such 
as these, where modulation is readily characterized by isolated peaks of enhanced and sup- 
pressed emission superimposed upon a quiescent state. We model the light curve with a 
sum of 10 Gaussians (one for each peak feature) and a constant (for the quiescent emission). 
This results in 31 free parameters. To obtain the best-fit parameter values, we perform 
chi-square minimization between the model and data. The final fit is shown as a solid line 
in Figure [1] and resulted in x 2 — 207.5 for 79 degrees of freedom. The large \ 2 is evidence 
that additional, small-scale variability is present in the light curve and/or that Gaussians do 
not ideally describe the flares. While the entire light curve is not ideally fit by this model, 
it provides a reasonable characterization of the variability and, more importantly, allows us 
to characterize the timing of the primary features in the light curve. 

We propose that the large peaks (MJD 54966.2, 55283.3, and 55601.3) and the dips 
(MJD 55303.4 and 55630.0) most accurately characterize periodic behavior in the light curve. 
The small precursor peaks that precede each large peak are unsuitable for measuring the 
period because they are observed to have variable shape, position (with respect to the large 
peaks), and height and because their significance above quiescence is small. The small mid- 
phase peaks are also unsuitable for measuring the period because the second peak, centered 
at MJD 55431.8, is poorly sampled, yielding an unreliable peak position estimate. Averaging 
the separations between peak centers, using large peaks and dips, we calculate a period of 
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321 days. We conservatively estimate the error to be ±5 days, the approximate average 
sampling rate during the large peaks and dips. Using this period, we have plotted the phase- 
folded light curve in Figure El This shows a strong similarity between modulation observed 
during different phases, particularly during the large peaks and dips. 

4.3. Significance of the Period 

To better estimate the significance of the 321 day periodicity, we performed Monte-Carlo 
simulations. By simulating many light curves with flare features similar in size and occur- 
rence frequency to those seen in the observed light curve, we can estimate the probability 
that the light curve is generated by a source that exhibits the observed flaring pattern every 
321 days, and not by a source that flares randomly in time. To characterize the flaring be- 
havior in the observed light curve, we again fit the 10 Gaussian + constant model to the light 
curve. Since we see multiple repetitions of each feature type (large, small, and mid-phase 
peaks and dips) in the light curve, it is possible to measure not only each feature type's rate 
of occurrence in the data set, but also estimate the range of values of each type's height 
and width by using the mean and variance. We define the rate of feature type occurrence 
as the number of occurrences per total time where the feature type may have, potentially, 
been observed. Given our observation cadence, the only time intervals where peaks or dips 
could not have been detected were during the date ranges MJD 54970 - 55050 and MJD 
55330 - 55424, where the source was too close to the Sun. During the entire period of Swift 
observations when peaks and dips could have been observed, we observed 3 large peaks, 3 
small precursor peaks, 2 dips, and 2 small mid-phase peaks. 

To generate the appropriate number of features of each type, we sample a random 
Poisson deviate corresponding to the calculated rate of occurrence. The shape of each feature 
is then created with a Gaussian peak of height and width generated by selecting a random 
deviate from a normal distribution with mean and variance measured in the corresponding 
feature type in the observed light curve. The mean position of each peak is generated by 
selecting a uniformly distributed random deviate on the observation baseline (MJD 54857 to 
MJD 55647). The peaks and dips are then combined with the best-fit value for the quiescent 
emission derived from data. The final simulated light curve is generated by sampling the 
model of superposed Gaussian peaks at the bin-center-times when Swift observed HESS 
J0632+057. We then calculate the \ 2 between the simulated light curve and the model 
that best fit the observed data. The fraction of simulated light curves with calculated x 2 
less than that of the data represents the probability of generating the observed light curve 
with a source that is stochastically flaring (with flares having the observed size and rate of 
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occurrence). In this way, we created 10 9 simulated light curves with randomly distributed 
flares and dips, none of which resulted in a \ 2 value less than that of the actual light curve, 
relative to the fit. Therefore, we find the false alarm probability for finding this periodic 
light curve from similar-sized stochastic flaring to be P < 1 x 10~ 9 . 
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Fig. 1. — The background-subtracted X-ray light curve of XMMU J063259. 3+054801 from 



the Swift-XBT observations in the 0.3-10 keV band. Figure 1(a) shows only the data and 
the time bin boundaries that were used for spectral fitting. Figure 1(b) shows the same data 
with the best fit of 10 Gaussians + constant model plotted as a solid line. Had the source 
been observable around MJD 55000, we hypothesize that it would have exhibited a dip state. 
To show this, we phase the best-fit model in the region of the MJD 55635 dip backwards 
by two periods and plot it as a dott ed line. We also ph ase forward the XMM data taken 
on 2007 September 17 (MJD 54360) (IHinton et al. 1120091 ) by two periods and plot its error 
(dominated by period uncertainty) as a gray shaded box. 



5. Discussion & Conclusions 



We have detected a 321 ± 5 day period in the 0.3-10 keV light curve of the unidentified 
TeV object HESS J0632+057. The observed periodicity has been shown to be significant 
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Fig. 2. — The ^-transformed discrete correlation function as a function of time lag, plotted 
with error bars (bottom) and divided by error bars (top), to give an estimate of significance. 
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Fig. 3. — The X-ray light curve of XMMU J063259. 3+054801 folded over the proposed period 
of 321 days. Zero phase has been arbitrarily defined as the date of first observation (MJD 
54857). The three phase cycles that result from this folding are designated with diamond, 
X, and square symbols, respectively. The lower panel shows the hardness ratio (2.0-10.0 
keV)/(0. 3-2.0 keV), folded over the same period and binned at 25 day intervals to improve 
the signal to noise ratio. The shown hardness data were fit with a constant (red dashed 
line), resulting in x 2 — 55.4 for 7 degrees of freedom, thus confirming variability. 
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by estimating the chance probability that the observed light curve periodicity is due to 
stochastic flaring of the source, resulting in a false alarm probability of P < 1 x 10 -9 . This 
implies a binary nature of HESS J0632+057, making it the fourth confirmed TeV binary 
(there are now 5 con firmed gamma-ray binaries, if one includes the recent announcement by 
Corbet et all toilh of GeV emission from 1FGL J1018.6-5856). The 321 ± 5 day period 
makes the system on e of the longest period Be star TeV/X-ray binaries and observations by 
Skilton et al. I (120091 ) show that HESS J0632+057 may have an unusually low luminosity. 



Figure [3] shows that the hardness does not vary significantly throughout most of the 
orbit, but it does reach a significant maximum when the light curve exhibits a dip feature. 
A constant line fit of the hardness data results in y 2 = 55.4 (with 7 dof), which shows that 
spectral variability is present, part icularly during the fl ux dip. This is consistent with the 



harder photon index measured by iHinton et al. I (120091 ) at a time corresponding to a dip 



feature in the light curve, phased forward by two periods (Figure [I]). This hardening of 
the X-ray spectrum during a flux decrease may be due to increased absorption of the soft 
X-rays in the source region, or it could be due to an orbital modulation of acceleration site 
parameters. If increased absorption due to a partial eclipse of the X-ray emission region 
is the origin of this spectral hardening, then it argues for a geometry of the binary system 
that allows the Be star and/or its equatorial disk to pass between the X-ray emission region 
and our line-of-sight. This would imply that the binary orbit is highly inclined and that 
the emission region is small enough to be eclipsed. However, this would likely be associated 
with an increase in N#, which is not observed in these data. An alternative, and more likely, 
explanation is that the spectral hardening during the flux dip is unrelated to absorption 
and instead caused by a change in acceleration site parameters, such as electron injection 
energies and cooling timescales, as a function of orbital phase. 

In addition to the flux dip feature discussed above, the dominant features of the periodic 
phased X-ray light curve are a large flux peak with a factor of 5-6 flux increase over quiescent 
flux and a moderate flux peak with a factor of ~ 2.5 flux increase over quiescent flux. These 
peaks are separated in time by about one half period. The flux dip immediately follows the 
large flux peak. Each of thes e featu res la sts roughly ~ 4 weeks. The recent increase in TeV 



flux reported by lOng et al. I (120111 ) and iMariotti et al. I ( 1201 ll ) coincides with the time of 



the recent large X-ray peak, which could imply related mechanisms such as synchrotron and 
inverse Compton emission. By analogy with known X-ray binaries, it is reasonable to assert 
that the peaks are the result of orbital modulation, but the system geometry and X-ray 
generation mechanism are not well understood. If the spectral hardening discussed above is 
due to absorption from the Be star and the surrounding region, then the geometry would be 
most easily solved if the large flux peak were due to periastron passage. However, studies to- 
date have not found evidence for optical radial velocity shifts in the optical counterpart MWC 
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148 (lAragona et al. 1120101 ; ICasares et al. 1 1201 if ). Additional, radial velocity measurements 



are required to understand the system's geometry 

While the light curve is not sampled identically in each orbital cycle, making orbit-to- 
orbit variability difficult to assess, it appears that some variation exists, particularly in the 
regions preceding each large X-ray flare. This may indicate the presence of other variability 
timescales relat ed to inhomogenei t y in the Be sta r disk/wind. Other X-ray binaries, e.g. 



LS 1+61° 303 (ISmith et al. 1 120091 ; iLi et al. 1 120111 ). show sporadic flaring behavior with a 



variety of timescales superposed on the binary orbit timescale, and this should be considered 
as a possible source of orbit-to-orbit variations for HESS J0632+057. Short timescale flaring 
analysis provides a way to probe the size of the emission region, as well as the power of the 
engine that must be feeding the associated acceleration site. If the dips indeed represent 
partial eclipses of the compact object by the donor star, this effect may also be used to 
constrain the size of the emission region and/or the size of a region of X-ray absorption. 

This work is supported at Pennsylvania State University by NASA grant NNX10AK92G, 
contract NAS5-00136, and grant NNX09AU07G. 
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